Physical Origin of the Boson Peak Deduced from a Two-Order-Parameter Model of 
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We propose that the boson peak originates from the (quasi-) localized vibrational modes as- 
sociated with long-lived locally favored structures, which are intrinsic to a liquid state and are 
randomly distributed in a sea of normal-liquid structures. This tells us that the number density 
of locally favored structures is an important physical factor determining the intensity of the boson 
peak. In our two-order-parameter model of the liquid-glass transition, the locally favored structures 
act as impurities disturbing crystallization and thus lead to vitrification. This naturally explains 
the dependence of the intensity of the boson peak on temperature, pressure, and fragility, and also 
the close correlation between the boson peak and the first sharp diffraction peak (or prepeak). 
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The dynamics of glass-forming liquids has been exten- 
sively studied in order to reveal the physical mechanism 
of the liquid-glass transition. One of the most myste- 
rious dynamic modes, which universally exists in many 
different types of glasses, is the low-frequency vibrational 
motion giving rise to the so-called "boson peak" in the 
inelastic scattering of light or neutrons 0. The boson 
peak originates from low-energy excitations having the 
energy of approximately 0.2-2 THz = 10-100 K and re- 
flects the excess density of vibrational states (DOS) in 
comparison with the Debye value. This excess DOS is, 
thus, also the origin of the low-temperature anomaly of 
specific heat Cp, namely, a hump in Cp/T 3 , which is 
again universally observed for different kinds of glasses. 
The physical origin of the low-energy excitations, or the 
low-temperature anomalies, remains to be clarified in the 
physics of condensed matter. 

The anharmonic contribution below 2 K was, on the 
other hand, successfully explained by two-level-systems 
models at least on a phenomenological level |^||, al- 
though the microscopic origin of a two-level-systems fea- 
ture and its universality are still not clear. The connec- 
tion of this behavior to the excess vibrational states is 
suggested by the soft potential model, based on the con- 
cept that both originate from the soft anharmonic nature 
of the potentials . Accordingly, the model predicts a 
close relationship between them. However, recent stud- 
ies by Sokolov et al. indicate the absence of such a 
correlation. A number of other theoretical models have 
also been proposed to explain the anomalous low-energy 
excitations commonly observed in glasses; e.g., the local- 
ization of vibrations due to homogeneous disorder whose 
spatial correlation extends to a medium range HI and 
the frequency resonance or localization of phonons due to 
nano-inhomogeneities (clusters or blobs) ]To|-p^|. How- 
ever, there has so far been no consensus on the physical 
origin of the boson peak. 

The fundamental questions associated with the boson 
peak, or the low-temperature anomalies, are as follows: 
(i) Why is the boson peak frequency several times 



smaller than the Debye frequency cjd of the system? 
What is the origin of such soft excitations? (ii) Is the 
anomaly produced by disorder effects? (iii) Is it related 
to the glass transition itself? (iv) How universal is it to 
glassy disorder systems? (v) Is it specific only to a su- 
percooled and glassy state or common to a liquid state? 
To answer all these questions, we need to clarify the true 
physical origin of a boson peak. One of the difficulties in 
specifying the mechanism of a boson peak is that the re- 
lationship between a boson peak and a liquid-glass tran- 
sition itself is not clear. 

In this Letter, we present a natural and coherent phys- 
ical explanation for the boson peak and its connection to 
glass-transition phenomena in the light of our two-order- 
parameter model of a liquid p?j[-p3t . We place a special 
focus on the experimental finding that a boson peak can 
exist even above the melting point T m for many strong 
glass formers such as Si02, B2O3, ZnCl2, and glycerol 
p4| p6),p8| . This means that a boson peak can exist even 
in an equilibrium liquid state for some strong liquids. 

First we briefly summarize our two-order-parameter 
model of a liquid |]l6|-^3|]. Our model is based on the 
concept that there are generally two types of local struc- 
tures in a liquid, normal-liquid structures and locally fa- 
vored structures, reflecting complex many-body interac- 
tions. The former is favored by density order parameter 
p, which attempts to maximize the local density and leads 
to crystallization, while the latter is favored by bond or- 
der parameter S, which originates from the symmetry- 
selective parts of the interactions and attempts to max- 
imize the quality of bonds. Locally favored structures 
are icosahedra for metallic glass formers, and tetrahe- 
dra and their ring organization for covalent-bonding glass 
formers. The locally favoured structure may have the lo- 
cal symmetry of a crystal, which is different from that 
of the equilibrium crystal, for a system with polymor- 
phism. Thus, we view a liquid as follows: locally favored 
structures with medium-range order are randomly dis- 
tributed in a sea of normal-liquid structures [see Fig. 1 
(a)]. It should be noted that a locally favored structure 
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has its own specific packing symmetry. A simple thermo- 
dynamic model of the two states with different degenera- 
cies [16 23 tells us that the average fraction of locally 



favored structures, S, is given by 
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exp[/3(AE - PAv)}, 
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where f3 — l/koT (fee: Boltzmann's constant), AE and 
Av are the energy gain and the volume increase upon 
the formation of a locally favored structure, respectively, 
and P is the pressure. Here g$ and g p are the degrees 
of degeneracy of the states of locally favored structures 
and normal- liquid structures, respectively. We assume 
that gs <C g p . This is the direct consequence of the 
uniqueness of locally favored structures and the exis- 
tence of many possible configurations of normal-liquid 
structures. Thus, there is a large loss of entropy upon 
the formation of a locally favored structure, which is 
given by As = h^\n{gs / 'g p ). The validity of eq. (|l|) 
was confirmed by the successful physical description of 
water's anomalies in terms of our model ]2l|-[23]|. The 
lifetime of a locally favored structure can be estimated 
as tlfs = T o exp(/3AG), where tq is the inverse of the 
attempt angular frequency and AG is the energy barrier 
to be overcome upon the transformation from a locally 
favored structure to a normal-liquid structure. 

Since the symmetry of locally favored structures is usu- 
ally not consistent with that of the equilibrium crystal, 
locally favored structures act as impurities for density or- 
dering, namely, crystallization. Thus, we can apply the 
knowledge of random-spin systems to this problem. Due 
to the random disorder effects of locally favored struc- 
tures, a liquid enters into the Griffiths-phase-like frus- 
trated state below T^, which is the melting point of the 
corresponding pure system and should be located near 
the real melting point T rn of the material (note that T m 
is usually located near the so-called mode-coupling T c 
^6|-|l9|). Namely, below the free energy of a system 
starts to exhibit a complex multi-valley structure. This 
leads to the non-Arrhenius behavior of the structural 
relaxation. The Vogel-Fulcher temperature To is then 
defined as a transition temperature from the Griffiths- 
phase-like state to the spin-glass-like nonergodic state. 

This physical explanation provides a satisfactory an- 
swer to the question of which physical factor governs 
the fragility of liquid |l(]-[ll|: In our model, fragility is 
determined by the strength of random disorder effects, 
namely, S. Thus, a liquid with larger S, which suffers 
from stronger disorder effects, should be "stronger" , or 
"less fragile" . This leads to a larger distance between 
and the Vogel-Fulcher temperature To for a stronger liq- 
uid, which is consistent with experimental results p6|-^9[. 
It is worth noting that the only way for a liquid to elimi- 
nate these random disorder effects, is "global crystalliza- 
tion" . While a system is in a liquid state, it is inevitably 
under the influence of frustration effects. 




FIG. 1. (a) Schematic snapshot image of a liquid. The 
black particles connected by bonds are participating locally 
favored structures, while all the other particles belong to 
normal-liquid structures. Note that the structure is rapidly 
fluctuating with time, (b) T-dependence of the boson peak 
intensity above T g , which is obtained by subtracting the 
quasi-elastic component from the total scattering. The raw 
boson peak intensity should be proportional to u)[n(u) + 1]S. 
Thus, the function I = Iq exp(AE/kBT) was fitted to the 
boson peak intensity scaled by ui[n(ui) + 1] by using 7o and 
AE/ks as adjustable parameters. Then ln(7/7o) is plotted 
against AE/k B T. AE/k B =3981 K, 3405 K, 1238 K, and 774 
K, for AS2O3 (open circles), Nao.5Lio.5PO3 (filled squares), 
B2O3 (filled circles), and glycerol (open squares), respec- 
tively. The data are taken from Refs. 24-27. The solid line 
is the theoretical prediction [see eq. (1)]. Note that below 
T g , S should become constant. Thus, the boson peak inten- 
sity scaled by the temperature bose factor becomes constant 
below T g in our model, in agreement with the well-known 
experimental finding. 

On the basis of this model, we seek the physical ori- 
gin of the boson peak. We assign the boson peak to the 
(quasi-) localized cooperative vibrational modes associ- 
ated with a locally favored structure and/or a cluster 
made of a few locally favored structures. Locally favored 
structures are characterized by a particular spatial sym- 
metry which is not consistent with the symmetry of the 
"equilibrium" crystal, and are distributed randomly in 
space. Their lifetime tlfs is much longer than that of 
normal-liquid structures, which are in a disordered state. 
Thus, a locally favored structure can have the coopera- 
tive vibrational eigen mode with an angular frequency of 
Wb P unique to that structural unit if tlfs is longer than 
l/wbp- This condition determines the high-temperature 
stability limit of the boson peak, Tb p . o->b P can roughly 
be estimated as <Jb P ~ 27ri> t /i?, where vt is the velocity 
of the transverse sound and R is the average character- 
istic length of a locally favored structure and a cluster 
made of a few locally favored structures. It is worth 
mentioning here that Uchino and Yoko Jl4|,|l5| specified 
the modes characteristic of particular types of medium- 
range ordering regions, which should correspond to our 
locally favored structures, for B2O3 and Si02 by ab initio 
molecular orbital calculations. 

In the light of this explanation, we now consider the 
questions (i)-(v) (see the introduction). The answers are 
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as follows: (i) A locally favored structure with medium- 
range order and its cluster are responsible for the boson 
peak and the low-temperature anomaly. Their charac- 
teristic size, R, is larger than that of the intermolecular 
bond, a. Thus Wb P is several times smaller than wd. (ii) 
Locally favored structures and clusters are isolated lo- 
cally ordered structures [see Fig. 1 (a)] and are not an 
average structure of supercooled liquid or glass, which 
is usually described by the term "medium-range order" . 
Thus, it may not be appropriate to simply say that the 
boson peak is due to disorder effects, (iii) It is a locally 
favored structure that causes frustration disturbing crys- 
tallization and is thus responsible for vitrification. In 
our model, thus, the boson peak is crucially related to 
the glass-transition phenomena, (iv) The boson peak 
should universally exist in any glass-forming liquid in 
principle since in our two-order-parameter model of the 
liquid-glass transition |jl6| |l9| it is the existence of locally 
favored structures that is the primary cause of vitrifica- 
tion, (v) Locally favored structures, which are responsi- 
ble for the boson peak, exist not only in a supercooled or 
glassy state, but also in a liquid state if Tb p > T m . 

Here we summarize important experimental findings, 
which any physical model of the boson peak has to ex- 
plain and thus can be used as criteria for determining 
the validity of a model, (a) The boson peak and the 
low-temperature anomaly of Cp are more pronounced in 
stronger glass formers ]7|j29| pl[|. (b) The boson peak in- 
tensity ]32| becomes weaker with increasing either tem- 
perature or pressure |3^,Q. (c) A boson peak has its 
own unique instability point, Tb p , which accompanies the 
appearance of quasi-elastic contributions. It is often lo- 
cated above the melting point T m for strong glass formers 

H-|||8), aS confirmed in Si °2' B 2°3, ZnCl 2 , and glyc- 
erol. These results seem to be contrary to the conjecture 
that Tbp ~ T c [ p7|j35| , although we cannot completely 
deny the possibility that they are exceptional cases. More 
importantly, they indicate that the boson peak can ex- 
ist even in an equilibrium liquid state above T m (i.e., 
Tbp > T m ). This cannot be explained by conventional 
models since they assume that an equilibrium liquid is 
completely homogeneous and has no specific structures, 
while a supercooled liquid has some medium-range order 
or specific disorder which is considered to be responsible 
for the boson peak, (d) There is much evidence [|36|— |3S|l 
which suggests the correlation between the boson peak 
and the first sharp diffraction peak (FSDP) (or prepeak) 
flof . Although there are experimental resultsjElJ con- 
trary to this prediction about the correlation ]36|, it may 
be true that the boson peak and the FSDP share some 
common origin in many systems. This is also supported 
by the fact that the FSDP has many features similar to 
the boson peak; for example, it often survives above T m , 
and its intensity decreases with increasing either temper- 
ature or pressure f42fl . 

An explanation of all the above-described findings (a)- 



(e) in a coherent manner at least on a qualitative level 
should be the first step toward an understanding of the 
physical mechanism of the boson peak. Thus, we at- 
tempt to define a coherent physical scenario rather than 
to quantitatively explain a particular experimental find- 
ing. Hereafter, we evaluate the validity of our explana- 
tion in terms of the above criteria (a)-(e). (a) According 
to our model, the magnitude of excess DOS, which is 
correlated with the intensity of a boson peak and the 
low-temperature anomaly of Cp, should be proportional 
to the number density of the locally favored structure, 
S. The two-order-parameter model of liquid |l6|-|l9|] in- 
dicates that stronger liquids suffer from stronger disor- 
der effects, which are characterized by larger S. Thus, 
the criterion (a) is naturally satisfied in our model. In 
our model |lq-|2C||, locally favored structures can be re- 
garded as impurities disturbing crystallization and this 
view is consistent with a recent conjecture by Quitmann 
and Soltwisch ]2S|] that the medium-range order respon- 
sible for a boson peak disturbs crystallization, (b) Our 
model indicates that the intensity of a boson peak |32] | 
may be constant below T g , but should decrease in pro- 
portion to the fraction of locally favored structures, S, 
in a liquid state above T g , with increasing temperature. 
Figure |l|(b) shows the T-dependence of the boson peak 
intensity for four strong glass formers, which is reason- 
ably explained by our prediction. The boson peak in- 
tensity should also decrease with increasing pressure in 
proportion to S above T g . However, there have been few 
experiments on the pressure effects in a liquid state. On 
the other hand, there are many experiments in a glassy 
state below T g , which clearly demonstrate the decrease 
of the boson peak intensity upon densification (see, e.g., 
Refs. 33 and 34). Such behavior is consistent with our 
model, but it may not be described by the thermody- 
namic relation, eq. ([l]), in an exact sense because of the 
intrinsically nonequilibrium nature of the densification 
process. Detailed studies on the temperature and pres- 
sure dependencies of the boson peak intensity "above T g " 
are highly desirable to more unambiguously confirm our 
prediction that the boson peak intensity is proportional 
to S between Tb P and T g . We could determine AE, As, 
and Av from such analyses, (c) In our model, a boson 
peak becomes unstable for WbpTLFS < !■ The tempera- 
ture of the instability of a boson peak, Tb p , is determined 
as 

Tbp ~ — AG/[fce ln(wbpTo)]. (2) 

Above Tbp, where tlfs < 1/^bp, the vibrational modes 
characteristic of locally favored structures should be over- 
damped. It is natural to expect a correlation between 
AE (the strong nature of a liquid) and AG. Thus, a 
strong glass former with large AG should have a high 
Tbp (sometimes above T m ). This is consistent with the 
observed behaviors p4]- |26| , p8[ | . Our model can also nat- 
urally explain the existence of the boson peak in a liquid 



3 



state above T m . A fragile glass former with small AG, 
on the other hand, should have a low Tb p . For a very 
fragile liquid, thus, the boson peak may not be observed 
in the liquid state. The fact that a boson peak survives 
above T m rules out the possibility that it is due to small 
clusters of the equilibrium crystal, (d) We propose that 
the FSDP or prepeak observed in the structure factor 
S(q) should reflect locally favored structures and their 
aggregates or clusters, in many cases. Recently, Mas- 
sobrio et al. [0 studied FSDP in GeSe2 by comparing 
their first-principles molecular dynamics simulations with 
neutron structure factors and found that the FSDP orig- 
inates from the medium-range order of tetrahedral sym- 
metry GeSe4, and its intensity decreases with increasing 
temperature. This is consistent with our explanation. 
Furthermore, Dzugutov et al. demonstrated that the 
FSDP is due to well-defined atomic-scale voids, which 
is associated with icosahedral local order for a simple 
monoatomic liquid. On noting that the locally favored 
structures for such a liquid are icosahedral structures, 
their finding is consistent with our explanation that the 
FSDP reflects the pseudo-periodic atomic density fluc- 
tuations associated with locally favored structures and 
their clusters. This explanation naturally leads to the 
conclusion that there should be a strong correlation be- 
tween the boson peak intensity and the intensity of the 
FSDP (or prepeak). Indeed, this is consistent with the 
observations [p6| |39| . Furthermore, this is consistent with 
the finding that in many cases the boson peak intensity 
and the intensity of the FSDP (or prepeak) |^2| both de- 
crease with increase of either T or P. However, it should 
be noted that the intensity of the FSDP (or prepeak) is a 
function not only of the number density of locally favored 
structures, S, but also their spatial arrangement. 

Finally, we provide other evidence to support our ex- 
planation. It is widely known that the addition of an- 
other component into a liquid changes its fragility. The 
famous example is a mixture of SiG"2 and Na20. Usually, 
Na20 is considered to act as a network modifier. We pro- 
pose, more specifically, that Na20 is the breaker of locally 
favored structures (probably, 6-member ring structures in 
the case of SiC^)- This concept naturally indicates that 
the addition of Na20 increases the fragility of liquids, 
and weakens the boson peak [45|| , the low-temperature 
anomaly of Cp G(|, and the FSDP, since it induces the 
reduction of the number density of locally favored struc- 
tures, S. The same scenario may also be applied to a 
mixture of B 2 3 and Li 2 @. 

We also mention a recent simulation study by Jund et 
al. p8[ , which systematically changes the strength of di- 
rectional interactions (AE in our terminology) for SiOv 
type liquids. It demonstrated that (i) an increase 
in the strength of directional interactions increases the 
intensity of a boson peak and the FSDP and also the 
number density of ring structures (locally favored struc- 
tures in our terminology) and (ii) a liquid with a more 



medium-range order is stronger, or less fragile. These 
findings are all highly consistent with our physical expla- 
nation. 

In summary, we propose a physical origin of the bo- 
son peak on the basis of the two-order-parameter model 
of liquid: Cooperative vibrational modes localized res- 
onantly to locally favored structures and their clusters 
may be the origins of the boson peak. This model can 
explain most of the experimental findings, at least on a 
qualitative level. We stress that locally favored structures 
exist not only in a supercooled or glassy state of liquid, 
but also in an equilibrium liquid state i/Tb p > T m . This 
is consistent with the experimental finding that the bo- 
son peak exists even above T m for some strong liquids, 
which is difficult to explain using previous models. 
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